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ABSTRACT 

A calibration is made for the correlation between the X-ray Variability Am- 
plitude (XVA) and black hole (BH) mass. The correlation for 21 reverberation- 
mapped Active Galactic Nuclei (AGNs) appears very tight, with an intrinsic 
dispersion of 0.20 dex. The intrinsic dispersion of 0.27 dex can be obtained if 
BH masses are estimated from the stellar velocity dispersions. We further test 
the uncertainties of mass estimates from XVAs for objects which have been ob- 
served multiple times with good enough data quality. The results show that the 
XVAs derived from multiple observations change by a factor of 3. This means 
that BH mass uncertainty from a single observation is slightly worse than either 
reverberation-mapping or stellar velocity dispersion measurements; however BH 
mass estimates with X-ray data only can be more accurate if the mean XVA value 
from more observations is used. With this calibrated relation, the BH mass and 
accretion rate may be determined for a large sample of AGNs with the planned 
International X-ray Observatory mission. 

Proper interpretation of the first AGN X-ray quasi-periodic oscillation (QPO), 
seen in the Seyfert galaxy RE J1034+396, depends on its BH mass, which is not 
currently known very well. Applying this relation, the BH mass of RE J1034+396 
is found to be 4+2 x 10 6 M . The high end of the mass range follows the rela- 
tionship between the 2f frequencies of high-frequency QPO and the BH masses 
derived from the Galactic X-ray binaries. 

We also calculate the high-frequency constant C = 2.37 M Q Hz~ x from 21 
reverberation-mapped AGNs. As suggested by Gierlihski et al., M B h = C/Cm, 
where Cm is the high-frequency variability derived from XVA. Given the similar 
shape of power-law dominated X-ray spectra in ultra-luminous X-ray sources 
(ULXs) and AGNs, this can be applied to BH mass estimates of ULXs. We 
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discuss the observed QPO frequencies and BH mass estimates in the ULXs M82 
X-l and NGC 5408 X-l and favor ULXs as intermediate mass BH systems. 

Subject headings: accretion, accretion discs — galaxies: individual: RE J1034+396 — 
X-rays: galaxies — X-rays: binaries 
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1. Introduction 

Active galactic nuclei (AGNs) are scaled-up versions of the Galactic X-ray binaries 
(XRBs), as recently evidenced by the physical link between the characteristic timescales 
and black hole (BH) masses corrected by the accretion rates (McHardy et al. 2006). The 
characteristic timescales are denoted by the break frequencies in the power spectral densities 
(PSDs). This suggests a similar broken power- law shape of PSDs from XRBs to AGN (e.g., 
Edelson & Nandra 1999; Uttley, McHardy & Papadakis 2002; Markowitz et al. 2003). 
The power spectra of XRBs also show multiple types of quasi-periodic oscillations (QPOs) 
other than the broadband noise. The first AGN QPO discovered in the Seyfert galaxy RE 
J1034+396 (Gierlihski et al. 2008a) stresses the similarities between the XRBs and AGNs, 
but its proper interpretation is hammered by its unknown BH mass. RE J1034+396 is 
an extreme Narrow-Line Seyfert 1 galaxy (NLS1), with the BH mass remaining largely 
uncertain (e.g., Bian & Zhao 2004). The BH mass determination is then a key issue since 
there are many controversies on the NLS1 mass estimates (e.g., Komossa & Xu 2007; 
Marconi et al. 2008; Zhu et al. 2009). 

We estimate the BH mass of RE J1034+396 from the X-ray variability amplitude 
(XVA), of ms (Nandra et al. 1997; Turner et al. 1999). A major difficulty in estimating BH 
masses from their XVA is the large scattering in the M B r— of ms relation (Lu & Yu 2001; 
Bian & Zhao 2003; Papadakis 2004; O'Neill et al. 2005; Zhou et al. 2007; Miniutti et al. 
2009; Nikolajuk et al. 2009). Note that the previous works were based on the BH masses 
derived from heterogeneous methods. This may cause additional scattering in the relation, 
since different methods of BH mass estimates may have quite different systematic errors. 
We therefore attempt to make a calibration of Mbh~ ff r 2 ms correlation based on only the 
reverberation-mapped BH masses; we further check this calibration result with BH masses 
determined from the dispersion velocity measures of these AGNs' host galaxies. We also 
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derive the high frequency constant as defined by Gierlihski et al. (2008b). This is useful for 
the BH mass estimates of ultra-luminous X-ray sources (ULXs) in nearby galaxies, since 
both ULXs and AGNs have characteristic power-law dominated X-ray spectra. 

This work is organized as follows. In Section 2, we re-analysis the XMM- 
Newton observations of RE J1034+396 and extract the light curves. In Section 3, we make 
a calibration of the correlation between the BH masses and XVAs and derive the high 
frequency constant. In Section 4, we discuss the uncertainties of mass estimates from XVA. 
In Section 5 we estimate BH mass for RE J1034+396 and discuss the implication of the 
AGN QPO. We find RE J1034+396 follows the relationship between the high-frequency 
QPO frequencies and BH masses derived from XRBs. Throughout this work, we use the 
cosmological parameters of h =70 km s" 1 Mpc -1 , Q m = 0.27, and f2 A = 0.73. 



2. XMM-Newton observations of RE J1034+396 

XMM-Newton observed RE J1034+396 on 2002 May 01 for 15.7 ks and on 2007 May 
31 for 93 ks. We screen the 93 ks data with the SAS v7.0 software. The X-ray events 
corresponding to patterns — 4 for the EPIC PN data are selected. We extract the source 
light curves from a 45" circle of the detected source position for PN, with the background 
being taken from six source-free regions with the same size. The presence of background 
flares and data gaps in the final 7 ks are checked and removed. 

Figure 1 shows the extracted 0.3 — 10 keV PN light curve (128 s bin -1 ), which is 
similar in shape to that in Gierlihski et al. (2008a). The light curve shows a periodic 
oscillation since t = 25 ks, with the frequency / = 1/3730 s ^ 2.7 x 10~ 4 Hz. Gierlihski 
et al. (2008a) have adopted the method proposed by Vaughan (2005) to quantify the 
statistical significance of periodicities against an underlying continuum which involves 
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dividing the periodogram by the best-fitting power law and using the known distribution 
of the periodogram ordinates to estimate the likelihood of observing a given peak. This 
method can be used to test the significance of candidate periodicities superposed on a red 
noise spectrum which has a power-law shape. The result indicates that the QPO is at ~5.6 
a significance level. 



XVA is the variance of a light curve normalized by its mean squared after correcting 
for experimental noise. For a light curve segment with iV bins (Nandra et al. 1997; Turner 



where Xj and Oi are the counting rates and uncertainties in each bin. \i is the arithmetic 
mean of the counting rates. 

Following O'Neill et al. (2005), the errors of of ms , which depend on the measurement 
uncertainties and the stochastic nature of the source, can be expressed as, 



where N seg is the number of available light-curve segments, <Tf rac is a fractional standard 
deviation, a frac = 0.74 for logM BH > 6.54 and a bac = 0.48 for logM BH < 6.54. A boot (of ms ) 
is the bootstrap uncertainty which comes from the bootstrap simulation accounting for the 
measurement uncertainties. 

XVA is calculated from a self-consistent analysis of an AGN's X-ray observation with 
a duration of at least 40 ks. The 2 — 10 keV light curves are subdivided into segments with 
a duration of 39.936 ks (binned at 256 s). The mean XVA is used if there are multiple 



3. Calibrating M Bli — of ms relation 



et al. 1999), 




(1) 
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segments. We take the measured XVA listed in O'Neill et al. (2005) if available. 



3.1. Reverberation-mapped sample 



Homogeneous BH masses are required for studying and calibrating the MBH-ofms 
relation. BH masses based on the reverberation-mapping (RM) method for 35 AGN 
(Peterson et al. 2004) and the small Seyfert galaxy NGC 4395 (Peterson et al. 2005) are 
currently available. Here we use the reverberation-mapped BH masses for the calibration of 
the M BH — af ms relation. We list in Table 1 the 21 reverberation-mapped AGN with at least 
one segment of X-ray data longer than 40 ks. 

In Figure 2, we show the correlation between XVA and the RM BH masses. We 
assume that there is a linear correlation between the M B h and XVA, y = a + (3x, where 
y = \og(MBn/M Q ), x = log(of ms ), with the measurement errors of e Xi for Xi and e Vi for y,i. It 
was shown that the Nukers' estimate is an unbiased slope estimator for the linear regression 
(Tremaine et al. 2002). The Nukers' estimate is based on minimizing, 



We apply the Nukers' estimate to derive the parameter of a = 4.85 ± 0.20 and 
(3 = —1.05 ± 0.08 for the 21 RM AGNs. The minimum x 2 P er degree of freedom is 1.6, 
indicating that either the uncertainties in the BH masses are underestimated or there is an 
intrinsic dispersion for BH mass estimates. 

To account for the intrinsic dispersion in the Mbh— of ms relation, we replace e y i by 
( e yi + e o) 1 ^ 2 ' where eo represents the intrinsic dispersion; eo is adjusted so that the value 
of x 2 P er degree of freedom is unity. This procedure is preferable if the individual error 
estimates e yi are reliable (Tremaine et al. 2002). Adding an intrinsic dispersion of 0.18 




N 



(3) 



-8 - 



dex decreases the value of x 2 P er degree of freedom to unity, and gives the best-fit results, 
a = 4.97 ± 0.26, and (3 = -1.00 ± 0.10. Then, 

M BH = lO 4 - 97±o - 26 (a r 2 ms )- LOO±o - lo M . (4) 

Note that the errors of the RM BH mass for IC4329A and NGC 4593 are quite large 
(Peterson et al. 2004). We exclude these two objects to ensure reliable error estimates of 
BH masses. The intrinsic dispersion of the M BB -a^ ms relation is 0.20 dex for the rest of 19 
AGNs. We also plot the objects with the BH masses estimated from the empirical virial 
relation; these objects (open circles in Figure 2) generally follow the best-fit correlation 
well. 

3.2. Sample with stellar velocity dispersions 

The key result we obtained is that the intrinsic dispersion of the Mbh— f r 2 ms relation is 
quite small, no larger than the uncertainties of RM BH masses of 0.3 — 0.4 dex (Peterson et 
al. 2004). We further check this result by using the BH masses estimated from the stellar 
velocity dispersions, since AGN hold the same relation between the BH masses and the 
stellar velocity dispersions (hereafter M B h— c* relation) as the normal galaxies (Gebhardt 
et al. 2000; Ferrarese et al. 2001). In Table 2, we list 21 AGNs with the stellar velocity 
dispersion measurements; only four objects are different from that listed in Table 1. 

Figure 3 shows the relation of XVA against the BH mass. The BH mass is estimated 
from the relation M BH = 1-349 x 10 8 (cr*/200) 402 M o (Tremaine et al. 2002). However, there 
is evidence for large scattering at the low BH mass end in the M B h-0"* relation (Greene & 
Ho 2006). We therefore use the RM mass of NGC 4395 (Peterson et al. 2005), and estimate 
the Mbh of MCG-6-30-15 from the empirical virial relation using the H/3 line (Zhou & Wang 
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2005). It is difficult to determine the actual errors of mass estimates in Table 2. Tremaine 
et al. (2002) claimed that the stellar velocity dispersion shows a Gaussian measurement 
error of 0.06 dex. Assuming no intrinsic dispersion in the Mbh - relation, the error 
of mass estimates from the stellar velocity dispersion should be 4.02 x 0.06 « 0.24 dex. 
Therefore, we add the error of ±0.24 dex for the BH masses in Table 2. 

We then apply Nukers' estimate to derive the parameter of a = 4.91 ± 0.24 and 
ft = —1.04 ± 0.09. The minimum x 2 per degree of freedom is 1.8, indicating that either 
the uncertainties in Mbh are underestimated or there is an intrinsic dispersion for Mbh 
estimates. We add an intrinsic dispersion of 0.27 dex to decrease the value of \ 2 P er degree 
of freedom to unity, and give the best-fit results, a = 5.15 ± 0.29, and ft = —1.01 ± 0.11. 
Then, 

Mbh = lO 5 - 15±o - 2 V™sr LO1±(U1 M . (5) 

The results are similar to that from the RM sample. This is a natural result since the 
RM method has been calibrated in agreement with the M B h-o"* relation (Peterson et al. 
2004). If the errors of M B h used here are underestimated, the intrinsic dispersion of the 
Mbh— ofms relation should be smaller. Since it is difficult to determine the actual errors of 
M B h estimates used here, we tend to use the calibrated relation with the RM sample, i.e., 
Equation (4). 

3.3. Weak dependence on Eddington ratio 

of ms may also depend on the second variable, Eddington ratio, rh E (e.g., Papadakis 
2004). Following McHardy et al. (2006), we hypothesize that logM BH = A + B\oga^ ms + 
ClogL bol . The fit is good with A = -(0.23 ±0.04), B = -(0.87±0.12), and C = 0.12±0.03. 
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Using rriE = Lboi/^Edd, then af ms oc m^. 14 / M-qh, indicating that the BH accretion rates only 
weakly affect the mass estimates. 

We further explore the dependence of of ms on the X-ray luminosity (e.g., Barr & 
Mushotzky 1986; Liu & Zhang 2008), FWHM of the broad components of the optical H/3 
lines (e.g., Leighly 1999), and the X-ray photon index (e.g., O'Neill et al. 2005). Results of 
partial correlation analysis show that the correlation between of ms and Mbh are the main 
correlation after removing the effects of other dependence (Table 3). 



Gierlihski et al. (2008b) have argued that the XVA at high frequencies can be used as 
an estimator of a BH's mass. Their hypothesis is that there is a universal power spectral 
shape for BHs at high frequencies (e.g., Lawrence & Papadakis 1993; Hayashida et al. 1998) 
above the PSD break frequency, f b , as illustrated in Fig. 4. This universal form is roughly 
a power law with an index a = 2, then P(f) = Cm(J / fr) 2 , where f r is an arbitrary 
frequency chosen to be f r = 1 Hz. 

Then XVA is calculated between frequencies f\ and / 2 , with f b < fi < f 2 , 



The assumption is that Cm is inversely proportional to the BH mass, Mbh = C/Cm, 
where C is the high frequency constant. Gierlihski et al. (2008b) calculated the C 
values from the Galactic XRB Cyg X-l. Different C values can be obtained if the data 
corresponding to the different spectral state of Cyg X-l are used; C = 1.24 M & Hz -1 for 
the low hard state and C = 2.77 M Q Hz -1 for the soft state. 

We calculate C from the 21 AGNs with the RM masses. We find C = 2.37 M Hz" 1 , 



3.4. Amplitude of High-frequency Variability 
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close to the value for the soft state of Cyg X-l. This means that C is not sensitive to M B h 
since the BH mass crosses over 8 orders of magnitudes from XRBs to AGNs. This also 
agrees with the fact that the amplitude of PSDs only slightly changes from XRBs to AGNs 
(see Fig. 18 in McHardy et al. 2004). Thus,, C can be reasonably assumed to be constant 
for AGN mass estimates. 

We emphasize that the C value from AGNs can be also applied for the BH mass 
estimates of ULXs in nearby galaxies, because the X-ray spectral properties of ULXs 
are similar to that of AGNs. Both ULXs and AGNs show power-law dominated X-ray 
spectra, together with a soft X-ray excess. Also, both ULXs and AGNs lack the 
disk-emission-dominated state as frequently observed in the Galactic XRBs. 

3.5. fi test and origin of intrinsic dispersion 

In Figure 5 we show the dependence of XVA on fc. If fi <C / 2 , of ms « C M f?/ fi, where 
fi corresponds to the lowest frequency of data sampled. That means, XVA is inversely 
proportional to f\ used. We calculate XVAs from a segment of 40 ks XMM-Newton data 
with a different f\. XVAs generally decrease with the increasing /i, but only a few AGNs 
follow the slope of —1. The measured deviation from of ms oc 1/ fi may be due to either a 
red-noise component in the light curve (Vaughan et al. 2003), or generally the assumption 
of a = 2, f b < fi <C f 2 is not always right for the sample. It was also known that the PSD 
of NGC 4051 shows a break frequency of 8 x 10" 4 Hz (McHardy et al. 2004), which is 
higher than the f\ used (for this work, f\ = 1/39.936 ks = 2.5 x 10~ 5 Hz). However, PSDs 
of most of AGNs in the sample are still unknown. These may contribute to the intrinsic 
dispersion of the M Bil — of ms relation. 
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4. Errors of M B h estimates from XVA 

The intrinsic dispersion of the M B h— of ms relation derived above is at a level of 0.2 — 0.3 
dex, no larger than the uncertainties of BH masses used here. Can the mass estimates from 
XVA be more accurate than those estimates from the empirical virial relation or from the 
stellar velocity dispersion? In this section we discuss the errors of the mass estimates from 
XVA. 

X-ray PSDs of AGN are not stationary on the timescales probed by typical XMM- 
Newton or Chandra observations. The red-noise nature of X-ray variability of AGN has 
been established (e.g., Lawrence et al. 1987; McHardy & Czerny 1987). For the red- noise 
process, XVA depends on the shape of PSD, and the time resolution and duration of the 
light curve. XVA may change from one part of the light curve to the next, even when the 
variability is produced by a constant PSD. Thus, different XVA values can be obtained 
if different data are used. Therefore, we attempt to estimate the uncertainties of mass 
estimates by looking at objects which have been observed multiple times with good enough 
data quality to make independent Mbh estimates. 

Figure 6 shows the distribution of XVA for the NLS1 galaxy MCG-6-30-15 and 
1H0707-495. MCG-6-30-15 was observed using XMM-Newton starting on 2000 July 11 and 
ending on 2001 August 4, with a total exposure time of ~480 ks. 1H0707-495 was observed 
using XMM-Newton starting on 2000 October 21 and ending on 2008 February 6, with 
a total exposure time of ~800 ks. We analyze all these data and filter the light curves 
to remove any background flares. We then subdivide the light curves into the 39.936 ks 
segments binned at 256 s. The XVA derived from each data segment shows the variability 
by a factor of 3 (~0.5 dex). This is roughly consistent with the error estimates by Equation 
(2). 

Gierlihski et al. (2008b) have analyzed RXTE Proportional Counter Array observations 
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of seven Galactic XRBs. The high-frequency variability derived from XVA changes by 
a factor of 2 — 3, even when the data corresponding to different spectral state are used. 
Therefore the error of mass estimates from a single data segment is slightly larger than 
the RM method of 0.3 — 0.4 dex, despite of the small intrinsic dispersion of the M B h— ofms 
relation. However, assuming XVA randomly scatters around the true value for the PSD 
of the process (Vaughan et al. 2003), the mean XVA of many data segments may reduce 
the error and thus the mass estimates can be more accurate than with the RM or stellar 
velocity dispersion method. 

We would also like to point out that the RM AGNs are a subset of BHs which are 
skewed toward low mass, high Eddington fraction of AGNs. These are probably the bulk 
of the AGNs for which reliable XVA can be computed with the present instrumentation. It 
would be not safe to attempt to apply this method to the objects like M31* or Sgr A* with 
the same calibration, since they both have extremely low accretion rates. The XVA can 
be a function of the photon energy. Thus the choice of spectral bandpass affects the XVA 
estimates. Nevertheless we do not address this issue in the current work since all the results 
presented here are based on the 2 — 10 keV band. However, if an X-ray instrument with a 
softer bandpass were used, or if one wanted to apply this method to higher redshift AGNs 
where both time dilation and redshifting of the X-rays were relevant, it would be necessary 
to develop an understanding of what would change. This will be investigated in a future 
work. 
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5. BH mass and QPO frequency 

5.1. BH Mass of RE J1034+396 

The X-ray variability of RE J1034+396 increases with the photon energy (Middleton 
et al. 2009). We extract the 2 — 10 keV light curve for the XVA calculation. The mean 
XVA of the two segments of light curves of RE J1034+396 is found to be 0.0234 ± 0.0081, 
corresponding to a BH mass of At\ x 10 6 M (See Figure 2). 

BH masses in AGN are usually determined from two methods: RM method and the 
Mbh — 0"* relation. We apply the empirical virial relation based on the RM method firstly 
and then compare the result from the stellar velocity dispersion, in order to estimate the 
BH mass in RE J1034+396. The gas in the broad line region (BLR) are thought to be 
virialized. The virial BH mass can then be calculated from 

M BH = j R B lr/ 2 FWHM H/3 /G', (7) 

where FWHM H/ g is the FWHM of the whole H/3 . For an isotropic velocity distribution, 
as generally assumed, / = y/3/2. -Rblr is obtained from the RM experiment or the 
£5100— -Rblr relation calibrated with the RM method. The BH mass in RE J1034+396 
calculated from Equation (7) is 6.3 x 10 5 M , when FWHM H/ 3 = 701 km s" 1 is used (Bian 
& Zhao 2004). However, recently it has been found that for NLSls, the broad H/3 line 
can be usually decomposed with a double Gaussian model, i.e., a very broad gaussian 
component (VBGC) and an intermediate gaussian component (IMGC) (Mullaney & Ward 
2008; Zhu et al. 2009). It is suggested that the BH mass should be calculated with FWHM 
of the VBGC (Zhu et al. 2009). Then, the virial BH mass should be corrected as 

M BHb = i?vBLR/ 2 FWHM| /3 (FWHMb/FWHM H ^) 2 /G', (8) 
where FWHMb is FWHM^ of the VBGC and -Rvblr is the radius of the very broad line 
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region (VBLR). -Rblr, is used as -Rvblr, because RM usually measures the size of the 
VBLR (Zhu et al. 2009). Its H/3 line can be well decomposed into three components as 
shown in Mason et al (1996): a narrow line component, VBGC and IMGC. Bian & Huang 
(2009) have recently analyzed the SDSS spectrum of RE J1034+396, by properly modeling 
the spectrum of its host galaxy. They found that the FWHM of the VBGC is 1690 ± 296 
km s _1 . We use Equation (8) to correct the BH mass, M BH b = (1690/701) 2 x 6.3 x 10 5 
M = 3.7 x 10 6 Mq. Considering the possible intrinsic scattering in the virial mass 
determination, the mass range of RE J1034+396 should be 2 — 7 x 10 6 M . This mass is 
in excellent agreement with the result from XVA, also roughly consistent with the range of 
1 — 4 x 10 6 M Q , determined from the stellar velocity dispersion of the host galaxy of RE 
J1034+396 (Bian & Huang 2009). 



5.2. QPO in RE J1034+396 

Figure 7 shows the scale of accretion flows between the BH mass and the QPO 
frequency across six orders of magnitudes. For a typical high-frequency (HF) QPO, 
f = 931(M BH /M ) Hz (Remillard & McClintock 2006), where /„ is the fundamental 
frequency. This relation is derived from three Galactic XRBs, XTE J1550— 564, GRO 
J1655-40 and GRS 1915+105, which display a pair of HFQPOs with 3:2 frequency ratios. 
The frequencies of XRBs are plotted for the stronger QPO that represent 2 x f . If the 
QPO in RE J1034+396 is a typical high-frequency QPO with / = 2/ or / = 3/ , it would 
correspond to a BH mass of 6.9 x 10 6 M© or 1.0 x 10 7 M , respectively. The high end of 
mass range of RE J1034+396 favors this AGN QPO as a HFQPO with / = 2/ . Middleton 
& Done (2009) suggested that the QPO is the AGN analog of the 67 Hz QPO in GRS 
1915+105, and not the HFQPOs seen in 3:2 ratios. However, that one scales the QPO 
frequency for the (4 x 10 6 /14) mass ratio between RE J1043+396 and GRS 1915+105, one 
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gets 76 Hz — within ~ 10% of 67 Hz, and well within the error bars on the BH mass for RE 
J1043+396. 

HFQPOs with the harmonic pairs of frequencies in a 3:2 ratio were often suggested to 
arise from some type of resonance mechanism (e.g., Abramowicz & Kluzniak 2001; Wang 
et al. 2007), where the frequencies are correlated with the BH mass and spin. The BH 
spin only changes the resonance orbits of accreting matter by a factor of a few. These 
may be the underlined physical reasons for the precise scale of accretion flows between the 
HFQPO frequency and BH mass. Rezzolla et al. (2003) proposed a model which requires 
that the two frequencies being oscillation modes of the same torus, rather than requires a 
resonance between the two QPO frequencies. The advantage of the Rezzolla et al. model 
for explaining what is seen here is that the higher frequency QPO is not seen, and the 
strengths of the different modes can vary independently. An alternative explanation for not 
seeing the higher frequency QPO might be that the PSD has not been made in the 10 — 30 
keV range where that QPO is usually seen in the Galactic XRBs, so this is not conclusive 
evidence in favor of the Rezzolla model over the resonance models. The QPOs previously 
seen are stellar mass BHs, which form through stellar evolution and grow insubstantially. 
They might all have quite similar spins. AGN, which can grow by both mergers and 
accretion, could have a different range of BH spin. Hence even if the QPO is due to the 
same mechanism as what produces the 3:2 frequency ratio QPOs in XRBs, a very realistic 
possibility is that the BH in RE J1043+396 has a mass closer to the most likely mass of 
4 x 10 6 M Q from the XVA analysis, but has a lower spin than the stellar mass BHs. 

5.3. QPOs in ULXs 

There are QPOs observed in the ULX M82 X-l (Strohmayer & Mushotzky 2003; Feng 
& Kaaret 2007) and NGC 5408 X-l (Strohmayer et al. 2007), but the BH masses of these 
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two ULXs were estimated ranging from 25 to 10 5 M (e.g., Dewangan et al. 2006; Wu 
& Gu 2008, etc). The non-standard accretion disk scenario (e.g., Soria et al. 2007) can 
explain the low soft X-ray temperature of M82 X-l without invoking the intermediate mass 
BH. However, by assuming the M82 X-l similar to the low-luminosity hard state of Galactic 
XRBs, Yuan et al. (2007) derived a mass of 9 - 50 x 10 4 M for M82 X-l by fitting the 
Chandra X-ray data together with the constraints of radio and infrared upper limits using 
the advection-dominated accretion flow model. Lower BH mass (10 4 M ) can be obtained 
if M82 X-l corresponds to the high-luminosity hard state. Assuming the QPOs in M82 X-l 
and NGC 5408 X-l are low-frequency QPOs, Casella et al. (2008) placed the constraints 
on the BH masses of 100 — 1300 M . Obviously, other possible QPO identifications will 
yield to different mass estimates. From the correlation between the spectral slope and 
the QPO frequency (known to exist in Galactic XRBs, see e.g. Vignarca et al. 2003), 
Fiorito & Titarchuk (2004) derived a mass of the order of 10 3 M for M82 X-l, although 
no such correlation was found in ULXs. Recently, Strohmayer & Mushotzky (2009) found 
that the X-ray spectral and timing properties of NGC 5408 X-l are quite analogous to the 
Galactic XRBs in the "very high" or "SPL" state. They also derived the high-frequency 
variability, Cm (Gierlihski et al. 2008b) for mass estimates. If we adopt the high-frequency 
constant C = 2.37 M Hz -1 from this work, their C M values correspond to mass ranges of 
3200 - 4000 M and 8400 - 10800 M for the two observations. This is not far from the 
scale between the HFQPO frequencies and the BH masses. 

6. Conclusions 

Calibration of the correlation between the black hole mass and the XVA shows that 
this correlation is quite tight, with an intrinsic dispersion of around 0.2 dex. The intrinsic 
dispersion may be caused by the uncertain shape of power spectral densities and the 
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stochastic nature of the red-noise light curves. The independent mass estimates from the 
X-ray variability derived from multiple observations show the uncertainties by a factor of 3 
(~0.5 dex). This is slightly larger than the errors of reverberation mapping or the stellar 
velocity dispersion measurements. However, the BH mass estimates can be more accurate 
if using the mean XVA derived from multiple observations. BH accretion rates only weakly 
affect the mass estimates. We find the high frequency constant C = 2.37M Hz -1 from 
the reverberation- mapped AGN. Given the similar shape of power-law dominated X-ray 
spectra in Ultra-Luminous X-ray sources and AGN, this can be applied to mass estimates 
of ULXs in nearby galaxies. The BH mass of RE J1034+396 is found to be 4+1 x 10 6 M 
from the XVA. The high end of mass range of RE J1034+396 follows the relation between 
the high-frequency QPO frequencies and the BH masses derived from the Galactic X-ray 
binaries, suggesting the precise scaling of accretion flows among them. We also discuss the 
observed QPO frequencies and BH mass estimates in the ULX M82 X-l and NGC 5408 X-l 
and favor the intermediate mass BH systems. 

We also find that X-ray variability as the mass estimator is independent of the 
luminosity, the distance and the viewing angle of the source; where luminosity is calculated 
by assuming the isotropy of emission. This is not true in many cases. Jets and outflows are 
common features and thus the beaming effect may play a role. The viewing angle may also 
change the inferred luminosity (Zhang 2005) considerably, since the X-ray emission from an 
accretion disk is generally not isotropic. Nevertheless, XVA provides quite reliable BH mass 
measurements without knowing the "true" luminosity of the source. Conventionally the 
BH mass in an AGN is often estimated from the measures of its broad line's width which 
may depend on the viewing angle of the source, unless the broad line region is completely 
spherically symmetric. For example, a narrower line profile can be obtained if it is seen from 
the pole-on angle and the BLR is flattened with a substantial Keplerian motion component. 
The radiation pressure may also affect the scale of the BLRs in AGNs. Apparently, XVA is 
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independent of these effects. Therefore, the method of estimating BH's masses with XVAs 
may be a good alternative for NLS1 galaxies and ULXs, whose BH mass estimates have 
remained largely uncertain so far. Long and deep X-ray survey observations may allow 
the BH masses and thus accretion rates determined for a large sample of AGN, with, for 
example, the planned International X-ray Observatory. 
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Fig. 1. — The source light curve extracted from the XMM-Newton PN camera in the 0.3 — 10 
keV band binned at 128 s for RE J1034+396. It is similar in shape to that in Gierlihski et 
al. (2008a). It shows a periodic oscillation of T = 3730 s since t = 25 ks. 
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Fig. 2. — Correlation between the X-ray variability amplitude and the BH mass. The filled 
points denote the 21 objects with the BH masses measured from the reverberation-mapping 
method. The best-fit linear relation given by the Nukers' estimate for the 21 AGN is also 
shown in the plot. The intrinsic dispersion of this fit is 0.2 dex. Most of objects show 
negligible scattering compared with the linear fit. The star denotes RE J1034+396. We also 
plot the objects with the mass estimated from the empirical virial relation; these objects 
(open circles) generally well follow the linear fitting relation. 
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Fig. 3. — Correlation between the X-ray variability amplitude and the BH mass. The BH 
masses are mainly estimated from the stellar velocity dispersion measurements, with an 
assumed error of 0.24 dex. See the text for the details. The Nuker's estimates give the 
similar results to that in Fig. 2, with an intrinsic dispersion of 0.27 dex. 
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Fig. 4. — The assumed power-law shape of the power spectral densities (PSD) at high 
frequency. The two dot-dashed lines denote the low frequency f\ and the high frequency f 2 
sampled. For this work, /i = 1/39.936 ks = 2.5 x 10~ 5 Hz, and f 2 = 1/2 x 256 s = 2.0 x 1(T 3 
Hz. To ensure the estimation of the integral over the power law PSD, the frequency break 
should be lower than fi. However, NGC 4051 has a higher frequency break of 8 x 10 -4 Hz. 
This may cause the deviation of NGC 4051 in the Mbh — ofms correlation. 
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Fig. 5. — Dependence of X-ray variability amplitude on j\ from a segment of 40 ks XMM- 
Newton data. of ms generally decreases with the increasing fi, but only a few AGN follow 
the slope of —1 (dashed line). 
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Fig. 6.— Distribution of X-ray variability amplitude derived from multiple XMM- 
Newton observations for the NLS1 galaxy MCG-6-30-15 (solid line) and 1H0707-495 (dotted 
line). 
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Fig. 7. — Relationship between the quasi-periodic oscillation (QPO) frequency and the BH 
mass. Three X-ray binaries, XTE J1550-564, GRO J1655-40 and GRS 1915+105 display 
a pair of high-frequency QPOs with a 3:2 frequency ratio. The frequencies of X-ray binaries 
are plotted for the stronger QPO that represent 2 x f Q . The line denotes the relation, 
f(Hz) = 1862(M BH /M )- 1 (Remillard & McClintock 2006). The observed QPO frequency 
in the ULX M82 X-1 and NGC 5408 X-1 is also plotted. See the text for details. 
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Table 1. Sample of AGN with reverberation-based mass 



Name 


z 




logafms 


Num. Seg. 


logLboi 


logm E 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


Mrk 335 


0.02578 


14.2±3.7 


-2.51±0.41 


1 


44.81 


-0.44 


PG 0026+129 


0.14200 


393 ± 96 


-2.88±0.66 


3 


45.97 


-0.72 


Fairall 9 


0.04702 


255±56 


-3.46 ±0.70 


6 


45.20 


-1.31 


Mrk 590 


0.02638 


47.5±7.4 


-2.51 ±0.13 


2 


44.76 


-1.02 


3C 120 


0.03301 


55.5 22 5 


-2.79 ±0.18 


4 


45.12 


-0.72 


Ark 120 


0.03230 


150 ± 19 


-3.42±0.91 


2 


45.18 


-1.10 


Mrk 110 


0.03529 


25.1 ±6.1 


-2.74 ±0.18 


1 


44.67 


-0.83 


NGC 3227 


0.00386 


42.2 ±21.4 


-1.62±0.16 


4 


43.33 


-2.40 


NGC 3516 


0.00884 


42.7 ± 14.6 


-2.43±0.10 


18 


43.83 


-1.89 


NGC 3783 


0.00973 


29.8 ±5.4 


-2.41 ±0.13 


8 


44.21 


-1.36 


NGC 4051 


0.00234 


1.91 ±0.78 


-1.06±0.09 


6 


42.88 


-1.50 


NGC 4151 


0.00332 


13.3 ±4.6 


-2.55±0.07 


29 


43.83 


-1.39 


* i oi 1 i 1 A Q 

r\j lZll-|-14o 


n nonnn 

u.usuyu 


1 A R _l_ A A 

140 ± 44 


Q i nor 

— o.Uo ± U.oo 


Q 

o 


4o. (V 


rt K£! 
— U.OO 


NGC 4395 


0.00106 


0.36 ±0.11 


-0.95±0.10 


6 


41.42 


-2.24 


NGC 4593 


0.00900 


5 36+ 9 37 

°- JO -6.95 


-1.85±0.33 


1 


44.03 


-0.79 


IC 4329A 


0.01605 


Q Qn+ 17 - 88 

y - yu -11.88 


-3.63±0.47 


6 


44.27 


-0.83 


Mrk 279 


0.03045 


34.9 ±9.2 


-2.63±0.36 


1 


44.83 


-0.81 


NGC 5548 


0.01718 


49.4 ±7.7 


-3.03 ±0.14 


16 


44.46 


-1.47 


3C 390.3 


0.05610 


287 ±64 


-3.41 ±0.41 


2 


44.90 


-1.66 


Mrk 509 


0.03440 


143 ± 12 


-3.24 ±0.59 


2 


45.23 


-1.03 


NGC 7469 


0.01632 


12.2 ± 1.4 


-2.33 ±0.27 


2 


44.67 


-0.49 



(1) Object name. (2) Rcdshift. (3) Black hole mass in units of 10 6 M Q , taken from Peterson et al. (2004) 

and Peterson et al. (2005). (4) Log of the mean X-ray variability amplitude. (5) Number of usable 
light-curve segments (39.936 ks per segment). (6) Log of bolometric luminosity. (7) Log of the Eddington 

ratio, rn E = L hoi /L E dd- 
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Table 2. Sample of AGN with stellar velocity dispersions 



Name 


z 


log(M BH /M ) 


(7* 


Ref. 


logofms 


Num. Seg. 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


Fairall 9 


0.04702 


8.36 


228 


4 


-3.46 ±0.70 


6 


Mrk 590 


0.02638 


7.84 


169 


2 


-2.51 ±0.13 


2 


Mrk 1040 


0.01665 


7.64 


151 


2 


-1.92 ±0.40 


1 


3C 120 


0.03301 


7.76 


162 


1 


-2.79 ±0.18 


4 


Ark 120 


0.03230 


8.32 


224 


1 


-3.42 ±0.91 


2 


Mrk 79 


0.02219 


7.38 


130 


3 


-2.15 ±0.42 


1 


Mrk 110 


0.03529 


6.99 


87 


1 


-2.74 ±0.18 


1 


NGC 3227 


0.00386 


7.45 


135 


1 


-1.62±0.16 


4 


NGC 3516 


0.00884 


7.95 


181 


1 


-2.43±0.10 


18 


NGC 3783 


0.00973 


7.65 


152 


4 


-2.41 ±0.13 


8 


NGC 4051 


0.00234 


6.53 


80 


3 


-1.06±0.09 


6 


NGC 4151 


0.00332 


7.23 


119 


2 


-2.55±0.07 


29 


NGC 4395 


0.00106 


5.56 


30 


1 


-0.95±0.10 


6 


NGC 4593 


0.00900 


7.30 


124 


2 


-1.85±0.33 


1 


MCG-6-30-15 


0.00775 


6.30 


159 


4 


-1.38 ±0.03 


48 


IC 4329A 


0.01605 


8.34 


225 


4 


-3.63±0.47 


6 


NGC 5506 


0.00618 


7.94 


180 


4 


-1.97 ±0.23 


2 


Mrk 279 


0.03045 


8.11 


197 


1 


-2.63±0.36 


1 


NGC 5548 


0.01718 


7.94 


180 


3 


-3.03 ±0.14 


16 


3C 390.3 


0.05610 


8.64 


268 


1 


-3.41 ±0.41 


2 


NGC 7469 


0.01632 


7.41 


133 


1 


-2.33 ±0.27 


2 



(1) Object name. (2) Rcdshift. (3) Black hole mass, mainly estimated from the stellar velocity dispersion. 
See the text for the details. (4) stellar velocity dispersion, in unit of km s" 1 . (5) references for stellar 
velocity dispersions, 1. Greene & Ho (2006); 2. Nelson & Whittle (1995); 3. Ferrarese et al. (2001); 4. 
Oliva et al. (1999). (6) Log of the mean X-ray variability amplitude. (7) Number of usable light-curve 

segments (39.936 ks per segment). 
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Table 3. Partial correlation analysis 



^2-lOkeVj M B H 








r MoX 


r LoM 




-0.85 


-0.69 


0.73 


-0.70 


-0.19 


FWHMH/3 , M BH 


fMa 


TFcr 


TMF 








— 0.85 


-0.52 


0.33 


— 0.84 


— 0.48 


I^-lOkcV, Mbh 




"TVcr 




»~Mcr,r 






-0.85 


-0.20 









r is the correlation coefficient. M, L, F and T denotes 
M B h, ^2-iokoV, FWHMH/3 and r 2 _iok e v, respectively. 



